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Abstract 



, The hyperconfluent third-order supersymmetric quantum mechanics, in which all the fac- 

^ , ' torization energies tend to a common value, is analyzed. It will be shown that the final potential 

as well can be achieved by applying consecutively a confluent second-order and a first-order 
SUSY transformations, both with the same factorization energy. The technique will be apphed 
^ I to the free particle and the Coulomb potential. 

cn ' 

^ : PACS: 03.65.Ge, 03.65.Fd, 03.65.Ca 

m 

o : 1 Introduction 

Nowadays, Supersymmetric Quantum Mechanics (SUSY QM) has became the standard technique 
for generating solvable potentials from a given initial one Moreover, it represents a powerful 

^ ! tool for designing Hamiltonians with a fixed prescribed spectrum (see, e.g, [7-9]). In its higher- 

\ order version, in which the differential intertwining operators are of order greater that one, it 

is well known that several seed solutions of the initial stationary Schodinger equation with an 
appropriate behavior are in general required for calculating the new potential, the eigenfunctions 
of the corresponding Hamiltonian, etc. [[T0l - fT8l . If for some reason just one of those seeds is 
available, one is driven to what could be called as hyperconfluent higher-order SUSY QM, in which 
all the factorization energies which are involved tend to a common value. 

In the past several works dealing with the confluent second-order SUSY QM have been elabo- 
rated [fT9] - l22] . Up to our knowledge, however, the hyperconfluent third-order SUSY QM has not 
been addressed explicitly. Of course, several papers involving the third-order SUSY QM have been 
done, but they are centered mainly in the case when the factorization energies are all different (see 
e.g. [|23l - |29l and references therein). 

In this article we aim to fill the gap by studying in detail the hyperconfluent third-order SUSY 
QM. In order to achieve this, we have arranged the paper as follows. In the next Section we 
will briefly review the confluent second-order SUSY QM. In Section 3 we will analyze the direct 
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approach to the hyperconfluent third-order SUSY QM, while in Section 4 we will address the 
corresponding iterative method. Section 5 explores the requirements that the seed solution has to 
obey in order to produce non-singular transformations as well as the eigenfunctions of the SUSY 
generated Hamiltonians. In Section 6 we will illustrate our general treatment by means of two 
specific examples, the free particle and the Coulomb potential. Our conclusions will be presented 
in Section 7. 



2 Confluent second-order SUSY QM 

Let us consider a one-dimensional Schrodinger Hamiltonian 

H, = -£, + V,{x). (1) 

The domain of definition or the corresponding system is denoted as "D = [x/, x^]. Thus, depending 
on the problem we are dealing with, and the consequent identification of xi and Xr, this domain 
could be the full real line, the positive semi-axis or a finite interval. The eigenfunctions and eigen- 
values associated to the discrete part of the spectrum of Hq, denoted by il)n{x)^ n = 0, 1, . . ., 
satisfy the stationary Schrodinger equation 

Hoi^n = -< + Vo^n = Eni^n, (2) 

as well as the boundary conditions 

1pn{xi) = IpniXr) = 0. (3) 

From now on we are going to suppose that all the eigenfunctions and eigenvalues of Hq are 
known. In the general formulation of the second-order SUSY QM one looks for a new Hamiltonian 

H2 = + V^ix), (4) 

which is intertwined with Hq by a second-order operator B2 in the way 

= B+Ho, (5) 

where 

Bt = i.-V{x)£ + I{x). (6) 

By plugging these expressions in the intertwining relationship ©, decoupling the resulting system 
of equations and solving it we arrive at [5 ,6 |: 

V2 = Vo-2r]\ (7) 

7 = Y + + , (8) 

T) = {ln[W{uuU2)]y, (9) 

where Mi, U2 are two seed solutions of the initial stationary Schrodinger equation associated to the 
factorization energies ei , 62 (in general different) 

HoUi = -u'l + VoUi = eiUi, « = 1, 2, (10) 
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and W{ui, U2) = uiu'^ — u[u2 denotes their Wronskian. Note that the seeds Ui, U2 could obey or 
not the boundary conditions of equation (|3]). 

The confluent second-order SUSY QM arises now as a limit procedure of the previous formal- 
ism when ei — )■ €2 — )■ e [19]. Note that, if the potential V2 is going to be different from the initial 
one, then ui and U2 cannot be just chosen as two linearly independent solutions of equation (fTOl) . 
since then W{ui, U2) = constant and therefore V2 = Vq. In order to produce non-trivial results, 
the right choice is to take ui as a standard eigenfunction of Hq while U2 becomes a generalized 
eigenfunction of rank 2 of Hq, both associated to e, namely: 

(ifo-eK = <=(K)-eK, (11) 

{Ho-e)u2 = ui {Ho-efu2 = U2 = {Vq - e)u2 - ui, (12) 

i.e., we are employing a Jordan chain of length 2. Expressing this in matrix language ElUl, this 
specific choice of basis M2} means that, in the restriction to the two-dimensional subspace of 
functions belonging to Ker(_B^), the initial Hamiltonian Hq is represented by a matrix (Hq) having 
a non-trivial Jordan structure of standard type: 

(i/o)=(j I). (13) 

Note that, if ui is given, then it is possible to determine the general solution U2 to the second- 
order equation (fT2] ) and to calculate then explicitly W{ui,U2) [21]. An alternative (and shorter) 
procedure runs as follows. First of all it is straightforward to show that 

W'{ui,U2) = U1U2 — U2U1 = —ul, (14) 

where we have used equations (nTl[T2]) . This implies that: 

W{ui,U2)=Wq- ul{y)dy = w{x), (15) 



XQ 



where xq G {xi,Xr). 

Let us note that, in order that the new potential 

V2 = VQ-2{ln[W{ui,U2)]}" = VQ-2[ln{w)]" (16) 

has not additional singularities with respect to Vq, then w{x) must not have nodes in (x;, Xr). This 
can be achieved by choosing a Schrodinger seed solution ui such that EOl : 



lim til = 0, ^- = ^liy) < 00, or (17) 

Jxi 

pXr 

lim Ml = 0, v^= I u\{y)dy< 00. (18) 

^^^^ Jxo 

With this choice, it tums out that w{x) becomes nodeless either for wq < —v_ in the first case or 
for Wq > Z/+ in the second one. Moreover, departing from the normalized bound states ipnix) of 

(2) 

Hq the normalized ones ipn {x) of H2 can be built up in the way: 

€>i.) = (19) 
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In addition, there is an eigenfunction of H2 associated with e which becomes as well square in- 
tegrable (we are using here a notation for this state which is appropriate for the purpose of this 
paper): 

nf 'ix) oc (20) 

Note that the confluent algorithm has been used to create bound states above the ground state 
energy of Hq [|20l . This possibility of spectral manipulation typically was outside the goals of the 
standard first-order SUSY QM. Moreover, the use of just one eigenfunction of Hq in the confluent 
case is advantageous compared with the second-order SUSY QM with ei ^ €2, which requires the 
knowledge of two Schrodinger seed solutions. 



3 Hyperconfluent third-order SUSY QM: direct approach 

In turn, let us analyze the hyperconfluent third-order SUSY QM, for which the three factorization 
energies converge to the same e-value, namely — e, i = 1, 2, 3. Similarly as for the second-order 
case of Section 2, we are going to use here a Jordan chain of length 3 of generalized eigenfunctions 
{ui, U2, U3} such that ui, U2 obey equations (|1 11121) while U3 satisfies 



{Hq - e)u3 = U2 



(Ho 







e)U3 - U2. 



(21) 



Equations (II 1I12I21I) mean that in the three-dimensional subspace of functions belonging to Ker(i?;^ 
this choice of basis implies that the matrix representing to Hq has a non-trivial Jordan structure of 
standard typeQ 




(22) 



Now, the hyperconfluent third-order SUSY partner Hamiltonians Hq and are intertwined 
by the third-order operator in the way 



BtH 



0) 



(23) 



where Hq is given by equation ©, H^ has the standard Schrodinger form 

H^ = -£, + V^{x), (24) 
and V3 is expressed in terms of the initial potential and the three seeds Mi, M2, M3 in the way: 

F3 = ^ - 2{ln[iy(Mi,n2,n3)]r, (25) 

with W{ui, U2, u-^) denoting the Wronskian of ui, U2 and M3 (we will give the explicit expression 
for in the next Section). A straightforward calculation leads to: 



W{ui,U2,U^) 



Ui U2 U3 
[ U'2 

u'{ U2 w 



U\ U'r, Un 



u'lW{u2, Us) - u'^W{ui, U3) + u'^W{ui, U2). 



(26) 



'Note that the matrices (Hq) of equations ( fT3] l and ( |22] | are non-hermitian despite Hq is hermitian. 
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By using now equations (II 1I12I21I) it turns out that: 

W{Ui, U2, Us) = UiW{Ui, Us) - U2W{ui, U2). 



(27) 



Recall that W{ui,U2) = w{x) was calculated in a simple way in the previous Section; thus, a 
similar procedure to obtain W{ui, U3) can be followed, leading to: 



W{Ui,U3) =Wi- 

J XQ 



ui{y)My)dy- 



(28) 



Given ui, and consequently the w of equation ([T?] ). it remains just to express U2 in terms of them. 
Let us note first of all that 



w = W{ui, U2) = Ml 

Henceforth 

U2 = ui /3i + 
Thus, a straightforward calculation leads to 



2 / 

Ul 



w{y) 



,■0 My) 

■' w{y) 



dy 



ui{y)u2{y)dy = wqI3i - I3iw{x) - w{x , . ^ 
XQ Jxo '^lyy) 



dy + 



f 


'w{y)' 


'XO 


.My). 



dy. 



(29) 
(30) 

(31) 



By plugging equations (I15I28I31I) into equation (ITTI) we arrive at: 



W{Ui,U2,U3) =Ui<fo 



XQ 



w{y) 

.My). 



dy \ = uif, 



with 



w{y) 



My) 



dy, 



fix) = fo- 

Jxo 

and fo = uji — wo/^i- Finally, the potential of equation (l25l) becomes: 

Vsix) = Voix) - 2{lnK(x)]}" - 2{ln[/(a;)]}' 
where f{x) is given by equation (|33l) . 



(32) 



(33) 



(34) 



4 Hyperconfluent third-order SUSY QM: iterative approach 

We are going to apply now two consecutive SUSY transformations departing from the initial 
Hamiltonian Hq: a confluent second-order one for generating V2 from Vq, which employs the 
two generalized eigenfunctions ui, U2 associated to e satisfying equations (II 1I12|) of Section 2; 
then a first-order transformation in order to obtain V3 from V2, which makes use of the general 
solution of the stationary Schrodinger equation of H2 associated to e. 

As for the confluent second-order transformation, we saw at Section 2 that the new potential 
V2 is given by 

V2 = Vo- 2{\n[W{u,,U2)]}" = Vo- 2[\n{w)]", (35) 

where the Wronskian W{ui,U2) = w{x) of the two generalized eigenfunctions ui, U2 of Hq 
associated to e is given by equation (fTSi) . 
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Concerning the first-order transformation, it turns out that H2 and H3 are intertwined by a 
first-order operator in the way: 



where H2 and are given by equations dH) and (|24] ) respectively, and 
with u^"^^ being the general solution of the Schrodinger equation 



(36) 



(37) 



eu 



(2) 



(2) 

From the results of Section 2 it is known that one solution is given by = ui/w (see equa 



tion (I20I) and 11201 |2T|| ). The other linearly independent solution is found by asking that 

W{u^^\u^2^) = 1 = [Mf^]^[4^V«f^]'' which immediately leads to = uf^ J^^dy /[uf\y)]'^. 
Thus, the solution m^^^ we are looking for to implement the first-order transformation takes the 
form: 



u 



(2) 



(2) , (2) 
CiU\ + C2U2 



Ul 



Cl 



-C2— < 



w 



C2 



xo 



My). 



dy 



(38) 



Hence, the final potential V3 resulting from applying the first-order SUSY transformation to the 
Hamiltonian H2, when using the seed solution given in equation (|38l) . becomes: 



n = V2- 2{ln[n(2)]}" = Vo- 2{ln[u,]}" -2 {in 



Cl 
C2 



XQ 




(39) 



Note that the two hyperconfluent third-order SUSY partner potentials Vi{x) of Vo(x) given by 
equations (|34l ) and (|39l) are exactly the same if it is taken /q = — Ci/c2. 

We can give, finally, the explicit expression for the third-order operator intertwining the 
initial and final Hamiltonians Hq and H-^ (see equation (|23l)): 



^3^2, 



(40) 



where 5^ and are given by equations <^ and (|37T ) respectively. 



5 Non-singular transformations and bound states of i/3 

As can be seen from equation (|25]) . in order that the potential V3 has no additional singularities 
compared with those of Vq, the Wronskian W{ui, U2, M3) given in equation (|32l ) should not have 
nodes inside V. This implies that both functions ui and / in this factorized expression should be 
free of zeros in this domain, in particular the seed solution ui which automatically leads to the 
restriction e < Eq, where Eq is the ground state energy of Hq. Moreover, for the second factor 
f{x) of equation (l33l) to be nodeless, the function w/ui should vanish either to the left edge xi of 
V or to the right one Xr- Here we are going to discuss in detail just the first case; the second one 
can be addressed in a similar way. 
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Let us choose first of all a nonphysical Schrodinger seed solution ui without nodes in V, 
obeying equation (fTTl) for e < Eq. Moreover, it is supposed that ui satisfies as well equation (fTTl ). 
Since w/ui should vanish for x — )■ a;/, we must have 

lim w = Wq + = 0, (41) 
which implies that wq has to be taken as 

Wo = -z/_ = - / ui{y)dy. (42) 



2/ 

U 
Jxi 

Therefore: 

w{x) = - u\{y)dy. (43) 

Jxi 

With this specific choice of ui and w, for most of the typical quantum mechanical problems it turns 
out that: 

Ui{x) 

Hence, the domain of the parameter /o such that f{x) is nodeless in (x/, x^) is given by 



w [x] 

lim — — — = and lim 

X^Xl Ul[x) X^Xr 



oo. (44) 



fo<-a^ = - I ^^dy. (45) 



X, ul{y) 



Let us note that in this /o-domain the third-order intertwining operator of equation POl 
tisforr 
the way: 



transforms the normalized eigenfunctions ipn of Hq into normalized eigenfunctions ipn^ of in 



^'^\x) = ^Mi^. (46) 
ViEn - e)3 

Moreover, the eigenfunction tpe^'^ of H3 associated to the eigenvalue e (compare equation (|38T)). 

Ha^) oc oc (47) 

turns out to be square-integrable in V, which implies that the spectrum of becomes 

Sp(i73) = {e}USp(ifo). (48) 

Note that, when /o — )■ — (t_, the hyperconfluent third-order transformation remains non-singular 
but the eigenstate ipl is not longer square-integrable. Thus, in this limit the two Hamiltonians H^^ 
and Hq become isospectral. 

On the other hand, for ui{x) being chosen as the normalized ground state eigenfunction iI)q{x) 
of Hq associated to Eq, it turns out that the previous equations (|4 111471) remain valid, the only 
difference is that now 

Jxi 

Thus, for /o satisfying equation ( |45] ). it tums out that e = Eq E Sp(//3), which implies that 

Sp(i/3) = Sp{Ho), (49) 
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i.e., the transformation is again strictly isospectral. However, when /o — )■ — ct- the ipe of equation 
(l47l) is not longer normalizable, meaning that in this limit e = Eq ^ 8^(11^), namely, 

Sp(i/3) = MHo) - {Eo}. (50) 

In this case, through the hyperconfluent third-order SUSY transformation somehow we 'delete' 
the ground state energy of Hq for generating H3. 



6 Examples 

Let us apply next the previous formalism to two physically interesting examples, the free particle 
and the Coulomb potential. 



6.1 Free particle 

The general solution of the stationary Schrodinger equation (fTTl) for the free particle with a negative 
factorization energy e = — /c^, A; > (for which Vo{x) = 0) is given by: 

= Ae'=" + 5e-'=". (51) 

In order to apply our method, let us use a nonphysical seed solution ui{x) satisfying equation (flTI) 
for xi = —00, i.e., let us make in equation (|5T1) B = and A = 1 so that: 

ui{x) = e^^ (52) 



With this choice, the calculation of equation P3l) leads to: 



2k 



Moreover, the evaluation of equation (|33T ) with Xq = produces 

- 



1 p2kx 

/(^) = /o + ^77^- (54) 



Note that this function does not have nodes for 

Hence, it is convenient to reparametrize this domain in the way: 



= w ^^^^ 

where Xi G (—00, 00). Thus, it is straightforward to show that: 

^k{x+xi) 

f{x) = cosh[fc(x - xi)]. (56) 

Finally, by plugging equations (|52I56I) into equation (|34l) . the hyperconfluent third-order SUSY 
partner potential of the free particle turns out to be: 

V^{x) = -2k'^sech.^[k{x - Xi)]. (57) 

This is the well know Poschl-Teller potential with one bound state at e = _Eo = —k"^, which has 
been also derived through first-order SUSY (see e.g. [0, page 30) and confluent second-order 
SUSY techniques (201. 
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6.2 Coulomb potential 

Working in spherical coordinates, separating the angular ones 9, (p, and making h = e = m = 1, 
the three-dimensional stationary Schrodinger equation for the Coulomb potential — e^/r leads to a 
one-dimensional problem characterized by the effective potential 

V.„(.) ^ ^(Itil, (58, 

where 0<r<oo,£ = 0,l,... The discrete energy levels En of Hq, for a fixed value of i, take 
the form En = —l/{n + £+ 1)^, = 0, 1, 2, . . . In order to apply our method, let us employ here 
the normalized ground state eigenfunction, 

1 / 2r ^ 



Mi(r) = , - — - e"^, (59) 

(£ + i)v/(2£ + i)! V^ + iy 

associated to the eigenvalue Eq = —!/{£ + 1)^. Let us start by calculating the w{r) of equation 
(l43l) with ri = 0, which leads to: 

7(a,x) being an incomplete Gamma function. Using this result and the expression for ui{r) of 
equation (|59l) it tums out that: 

which vanishes for r — )■ 0, as required. The calculation of the /(r) of equation (|33l) with ro = 
produces now: 

, _ , >^ >^ 7(fc + m-2£-l,g,) 

fc=2^+3 m=2^+3 

7(2£ + 3,^) , / 2r 



2r(2£ + 4) - - \^ ' ' ' ' ' £ + 1^ 
(£ + 1)^ 7(m + 2£-F5,g,) ' 

4 ^ (m + 2)(m + 2£ + 3)!' ^ ' 

m=0 ^ ' ^ ' 

which is nodeless in (0, oo) for /o < 0. The hyperconfluent third-order SUSY partner of the 
effective potential (l58l) becomes finally: 



f{r)ul{r) 



(63) 



where ui{r), w{r) and /(r) are given by equations (I59I60I ) and (|62l ) respectively. 

The first two terms of equation (|63T ) correspond to an effective potential different from the initial 
one (compare equation (|58T)). This difference is also reflected in the energy levels of a potential 
composed only of these two terms, which are given by En = — l/(n + £+ l)^, n = 1,2,... Thus, 
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it is natural to interpret that the third term of equation (|63l ) is the main responsible of supporting 
the ground state energy of V3(r) at Eq = —!/(£+ 1)^. 

Let us note that the family of hyperconfluent third-order SUSY partner potentials given by 
equation (l63l) is different from the ones which have been derived either by first-order SUSY [115] 
[T6ll30l - [32l or by second-order SUSY transformations [Il5] [16] |211 (just compare the centrifugal 
terms of each family). 

In particular, for £ = it tums out that, departing from the Coulomb potential without centrifu- 
gal term, Vo(r) = — 2/r, we arrive at a new one-dimensional potential with a non-trivial centrifugal 
term given by 



V,(r) 



2 2 



f{r)ul{r) 



where now 



uAr) 



2re 



1 



wir] 



7(3, 2r) = (2r2 + 2r + l)e" 



-2r 



1 1 °° 

fir) = fo- Y^7(3, 2r) .F^ (1, 2; 3, 4; 2r) + - 



m=0 



7(m + 5, 2r) 
[m + 2) (m + 3)!' 



(64) 

(65) 
(66) 



As an illustration, the isospectral potentials Vo(r) = —2/r and the V3(r) of equations (I64II66I) for 
/o = —1/10 as functions of r are shown in figure 1. The corresponding energy levels of H^, and 
Ho are given by En = —l/{n + 1)^, n = 0, 1, 2, . . . 




Figure 1 : Coulomb potential Vq (r 
partner Vs^r) given by equations 



2/r (gray curve) and its hyperconfluent third-order SUSY 
for/o = -l/10. 



Let us note that the well arising in V3(r) at the neighborhood of r ^ 1 is induced by the 
third term of equation (|64] ). which is also the responsible of supporting the ground state energy at 
Eq = —1. Let us recall that this level was not present in the effective potential of equation ([58]) 
with£ = 1. 



7 Conclusions 



In this paper we have addressed the hyperconfluent third-order SUSY QM through two different 
(but equivalent) approaches, namely, direct and iterative one. It was found the explicit expression 
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for the Wronskian, the most relevant quantity which determines the form of the new potentials, the 
eigenfunctions of the associated Hamiltonians, etc. 

The requirements for the seed solution to produce non-singular SUSY transformations were 
as well explicitly determined. Note that, from considerations taking into account the order of the 
transformation, through the hyperconfluent third-order SUSY QM one obtains a three-parametric 
family of potentials (for a fixed factorization energy). However, since we had to impose two 
requirements on the solutions employed in the iterative approach, it turns out that the non-singular 
potentials for r G (0, oo) belong just to a one-parametric subset of the general three-parametric 
family which one is able to build up. 

Our general procedure was illustrated by means of the free particle and the Coulomb poten- 
tial. In particular, the last case illustrates clearly that the non-singular one-parametric family of 
potentials derived through the hyperconfluent third-order SUSY QM is different either from the 
set which can be achieved from a first-order SUSY transformation I115[|16[|30] - |321 or from the one 
which can be generated through the confluent second-order SUSY transformation [|2n . 
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